Globally controlled universal quantum computation with arbitrary subsystem dimension 
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We introduce a scheme to perform universal quantum computation in quantum cellular automata (QCA) 
fashion in arbitrary subsystem dimension (not necessarily finite). The scheme is developed over a one spatial 
dimension A'-element array, requiring only mirror symmetric logical encoding and global pulses. A mechanism 
using ancillary degrees of freedom for subsystem specific measurement is also presented. 
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Quantum computation involves control at some level of a 
large number of quantum subsystems. The traditional circuit 
based approach is to have a set of subsystems for storing quan- 
tum information, denoted q-sites, which are fully addressable 
not only individually but in subsets [1]. It is known that uni- 
versal quantum computation is possible, independent of sub- 
system dimension, when arbitrary single q-site unitaries and 
2-q-site entangling gates are available [2, i]. Yet addressing 
individual subsystems in large arrays is extremely challeng- 
ing and can impose significant errors due to miss-alignment 
of the fields which unintentionally act on neighbors to the tar- 
get (cross-talk) or miss the target. Thus as a mean to avoid it 
another option appeared: global control schemes [4, 5, 6]. 

The philosophy behind global control is to reduce the in- 
teraction with the array of q-sites to require only global ma- 
nipulation, implemented for instance through global fields ho- 
mogeneously coupled with all q-sites. In general such global 
control schemes require a natural evolution (time step) for the 
array, and tailored sequences of global pulses which translate 
physical asymmetries in the array into control of particular 
sites or, equivalently, chronological control into spatial con- 
trol. When the resultant evolution is a set of gates which act on 
small neighborhoods in parallel, this is also a quantum cellular 
automata (QCA) model. The models for global control have 
been so far concerned only with qubits [ '^], however with 
the development of higher dimensional computational models 
(using qudit computation and continuous variables(CV), also 
called qunat computation) that show advantages in terms of 
efficiency and robustness [ ], the natural direction is to find a 
way to implement such models in a globally controlled fash- 
ion. 

The aim of this paper is to develop such a model in one 
spatial dimension, inspired by a previous protocol restricted 
to qubits [6] and recent results on globally controlled trans- 
port of qudits and qunats [' .]. The main technical difficulty is 
the presence of more complex phases that appear in higher di- 
mensions, in contrast with the {1,-1} phase in the qubit case, 
and solving equations for discrete variables which are defined 
modulo the dimension d of the logical q-site under consider- 
ation. Fortunately, in [8] we developed most of the tools we 
will need as well as revising some of the known results avail- 



able in the literature [9]. 

The paper is organized as follows: in Sec. I we review some 
mathematical results on bases for the operator space particu- 
larly exploring a Hermitian basis so that we can obtain uni- 
versal QC through gates generated by physically accessible 
Hamiltonians. In Sec. II we develop our scheme and show how 
certain sequences of global pulses plus a natural time step, act- 
ing on a N-site array with a mirror-symmetric logical encod- 
ing, can generate single site arbitrary unitary gates as well as 
nearest neighbor entangling gates and thus achieve universal 
quantum computation in a quantum cellular automata fashion. 
Input and output of information into the array is discussed in 
Sec. Ill and we conclude with a summary of the results. 



I. SINGLE QUSITE QC 

The state spaces for discrete and continuous q-sites are 
qualitatively different as unitary operators on these spaces 
are generated by the algebras su(c/) and polynomials in 
{^,^,1} respectively where the canonical commutation rela- 
tion [q,p] = ih is satisfied. Motivated by results in Ref. [II] 
we can chose a set of generators for single q-site quantum 
gates using Hermitian counterparts of Weyl pairs which sat- 
isfy the the commutation relations of the generalized Pauli 
group for qudits and the Weyl representation of the Heisen- 
berg commutation relation for CVs. 



A. The qudit case {d finite) 

The state space for a single qudit is yi^i = ^P^^c{\j)}'jZo 
and since the global phase is irrelevant, unitary gates are ele- 
ments of the group SU [d). A unitary operator basis on this 
space is {X{a)Z{b);a,b e Z,/} where X = ^^^o l-^+l) (■^l- 
and Z = Tf'Zo ^'^^ 1*) with additional modulo d satisfying 
X'' = l =Z''. We adopt the notation X (a) =X" and Z(Z7)=Z^ 

The group commutator is Z{b)X{a) = e~3'"^X(a)Z{b) . How- 
ever if we are to perform quantum gates generated by Hamil- 
tonians we require a Hermitian basis of 5u{d) such as the fol- 
lowing. 
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B{a,b) = e'^'--t'X{a)Z{b)+e-'^"-bZ{-b)X(-a) 
= co^{<!^a.b){X{a)Z{b)+Z{-b)X{^a)) 
+ ?.m{(Sf,,b)i{X{a)Z{b)-Z{-b)X{-a)) 



(1) 
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where (|)„^ = K/4 + Kab/d. Since the two Hermitian sum- 
mands commute we can write an arbitrary unitary U £ SU (d) 

as 

U — lla,f,=o'^ /ON 

ga„.4sin(0„,,)(X(a)Z(fo)-Z(-fo)X(-fl)) 

where the aa.b are real numbers (note the presence of c/^ — 1 
independent real parameters). It follows that if we can per- 
form e™~ {x{a)z{b)±z{-b)x{-a)) ^ f^j. ^^j^^g ^Yien we are able 
to construct an arbitrary unitary. 

B. The CV case {d ^ oo) 

The state space for a single qunat is J-Ccv = Con- 
sider the following operators on this space: 

Z(a) = e'""? and X(p) e-'P'\ (3) 

satisfying Z(a)X(p) = e'"Px(p)Z(a) which is the Weyl rep- 
resentation of the Heisenberg group (having set Ti~ 1). We 
can fix a computational basis |^) satisfying the relations 
Z(a) 1^) = e'"'? 1^) and X(p) |^) = |^ + p). Unitaries gener- 
ated by the set of Hermitian operators corresponding to the 
real and imaginary parts of X (P)Z(a) can be constructed, and 
in Ref. [ 1 1 ] it is argued that such gates are universal for single 
CV computation. To elucidate this point, consider the unitary 
evolutions obtained using the gate library 

|g/«cosa),,3 g/6sina),,5 g/.costO;^ g,ysinco,,^|^ fl,Z,,C,t/ £ R (4) 

we shall argue that this components are enough to have sin- 
gle qunat universal computation i.e. are enough to generate 
any Hamiltonian evolution. The core of the argument is a re- 
sult from functional analysis [12] which shows that the set 
{sin«x,cos«x} is a basis for the functions space, in particular 
to the polynomials space. 

oo 

/(x) = Co + ^ C„cos«.x + D„sinnx, V/(x) (5) 
11=1 

It follows then that we can approximate any Hamiltonian in 
q and p through a convenient sequence of basis elements ac- 
tions, or products of sums, with arbitrary accuracy. However 
this ability comes at a cost, we are able to reproduce the func- 
tion only within a domain {— L,L}, and then it starts repeating 
itself effectively making the phase space periodic q = q + riL, 
p = p + ml! for some L,L' and all n,m, allowing us to define 
periodic coordinates qL,PL- 

q,p > qL,PL 

[qj]^i — > [qL,pL]=i 

In this way we would have, for instance 



within {— 7t,7t}. 

Since the terms inside the argument of the exponential com- 
mute we can apply then a product of exponentials and truncate 
at some «max to approximate a polynomial in q: 

"omax _ ^emax 
n=l n=\ 

where (n) is the discrete cosine(sine) Fourier transform 
of /(^l) [ ' "']■ Given the ability to generate an arbitrary qubic 
polynomial in qi and a quadric polynomial in pi, it is possi- 
ble using the Lie-Trotter theorem to generate any polynomial 
fiqiipL) [ ^]- This allows us to have then the same results 
in the standard phase space, namely Hamiltonians and their 
respective eigenstates written in terms of the new qL,PL co- 
ordinates. Alternatively one could directly span the any or- 
der polynomial as an expansion of sine and cosines. Thus, 
in principle, universal single qunat gates with arbitrary accu- 
racy are viable in the {— L,L} domain with currently physi- 
cally achievable Hamiltonians. Thus to approximate a general 

Hamiltonian we must be able to construct Hamiltonians of the 
form ^j- g/(x(«)z(i)-z(-i)x(-«))^ 

so we start from the elements in (4), now if we can also per- 
form global e'('^ +P rotations we can induce the transforma- 
tion 

q ^ q cos CO + /3 sin CO 
p — ^sinco + /3cosco 

and thus e-'(9^+P^)™e'"'=o'*{P4)e'('?^+P^)™ = gia-cos{aq+bp) ^-^^ 
a,b Ij, choosing cosco = a/y/a^ + b^, sin CO = b/^/a^ +b^, 
and p = VcP^+l?. Similarly, conjugating evolution gen- 
erated by sin(p^) we can obtain e'«sin(«9+6p) Hence 
using the fact that those evolutions commute and that 
X{a)Z{b) + Z{-b)X{-a) = 2{3i{e'"'') cos (ap + bq) - 
3 {e'"'') sin (ap + bq)) and X{a)Z{b) - Z{-b)X{-a) = 
2(3(e'"*)cos(fl;9 + ^7^) + 3i{e'"'') sin {ap + bq)) we can 
construct any unitary U = or 
jj _ ga(x(a)z(b)-z(-b)x(-a)) ^ p^j. universal qunat computation, 

the authors in Ref. [3] showed that it suffices to have arbitrary 
single qunat gates and at least linear coupling in q between 
qunats. Hence we confirm the argument in Ref. [11]. 

In summary, to implement a universal scheme for quan- 
tum computation, we will need to be able to apply the set 
of operators ( 1 ) on any site s, and furthermore we need an 
entangling gate between at least nearest neighbors {s,s + 1}. 
In the following section we shall explore a QCA scheme on 
a 1 — D lattice with N-sites, where using only global pulses 
we can achieve the required addressability for arbitrary di- 
mension (d) elements. Most of our calculations are equiva- 
lent as we are mainly using the algebraic and transformation 
properties of the Z{u)X{v) operators, thus a unified notation 

Z{u)X{v) = C'"'X(v)Z(m) where ^ e'^ for finite c/ and ^ = e' 
for CV, is in order. 
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II. UNIVERSAL COMPUTATION WITH GLOBAL 
CONTROL 



Inspired in the scheme developed by Raussendorf [6] we 
generaHze it to arbitrary dimension. Anticipating the results 
of the the paper we shall fix at this point a mirror-invariant 
encoding of the initial state, effectively turning Eq. (6) into 
a single logical site gate. We will encode M logical q-sites 
into N physical sites, the mirror symmetric encoding leads us 
to two somewhat similar encodings depending if is odd or 
even. If is odd we can encode M ~ logical q-sites, 
while if is even we can encode M = N/2 q-sites i.e. q- 
sites encode [A^/2J logical q-sites. More specifically we do 
this by mapping pi,...,M ^ Pi,...,m «> Xmiddie «> Pm,...,i, where 
the Xmiddie Contains no encoded information and is dropped if 
A^ is even. 

We start recalling the step operator T for arbitrary dimen- 
sion introduced in [8] for a A^-elements chain, 

N N-l 

where Fj is the Fourier gate acting on element j (satisfying 
F'^ = 1) and CZ, ,+i is the generalized control phase gate be- 
tween sites ; and ; + I. For qudits Fj = ■;^Lr.7=o^" I'') j 

andCZ,i+i = I^^'^Z' ®Z^' while for CV, F = e't('''+''') and 
CZ, ,+i = e'm®ii+\ ^ Henceforth, all operations are global so 
we write F = Yfj=\ Ff^ CZ = n£i' CZ,- ,+i . The compo- 
sition F^T^^^ is a reflection about the middle on the state of 
an A'-element chain [8]. Consider the operator. 



(a) 



PieN+i)TP{eN)TP{eN-i)...TP{eo) 
P{eN+i)iTP{eN)T-'){T^PiEN-i)T-^) 

...(r^+'p(eo)r-(^+>))r^+i 

N+l 
111=0 



where P{£s) = (Sl^li Xi{—es)Zi(Es), and also the following ac- 
tion 



y(a,«,V,Z) = (f'2f^)-lj-J^-/p/2(X(„)Z(v) + (X(„)Z(v))1-) 

j 

(/r2f^-)-j-Jg-'-p/2(X(»)Z(v) + (X(„)Z(v)) + ) 
j 

(F2fjF2-)-j-Jg;p/2(X(»)Z(v) + (X(»)Z(v))t) 

j 

= ■ J~Je'«/2(X(»)Z(v) + (A-(«)Z(v))t)A,-,^ 



(6) 



for some /, where we have defined A{j,l) = 8/,/ + 8j_N+i-i, 
the symmetric delta function. Here we have used the fact that 
for mirror symmetric operators (which is the only kind used 
in our control), 7^2^-^+! ^^^^ identity. The overall result 
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FIG. 1 : (a) The architecture consists of a chain of N physical q-sites 
encoding M logical q-sites (shown here for A' even) in a mirror sym- 
metric pattern, i.e. identifying site / with site A' + 1 — Z. To perform 
readout, an ancillary q-site (another qudit or another CV degree of 
freedom) is associated with each data q-site shown here as a another 
parallel chain, though the ancillary and data q-sites may be spatially 
co-located. The ancillary q-sites are initialized in |0) *''*'. (b) A circuit 
for executing the SWAP gate between data and ancillary q-sites at lo- 
cations / and A' + 1 — / . We require the ability to execute a global gate 
CZj g between the two chains, where we recall that the global swap 
gate can be executed through CZ and global F gates. Also note that 
the remaining sites, on which no F gate is performed, undergo trivial 
evolution. Finally, the required sites are swapped and we can per- 
form a global measurement on the ancillary chain to readout logical 
q-site /. The process can be repeated subsequently for the remaining 
sites. Note that only global operations/measurements are used. 



is then a unitary generated by X{u)Z{v) + (X(m)Z(v))^ on site 
I and its mirror image only. Here where u,v G {0,...,d— 1} 
for the qudit case and e M for the CV case. Before proceeding 
with the calculation it is worthwhile giving some motivation to 
such claim: if one does the calculation for some short chain, 
say 6 qudits (or qunats), and sees how a homogenous P(e) 
evolves using the time step introduced in [;s] one realizes that 
every site (and the mirrored position) ends up with a differ- 
ent element of the Clifford group, and thus a Clifford operator 
commuting with such time-evolved operator will gain a differ- 
ent phase which will depend on the particular site of the chain 
it is acting. This was first exploited for qubits in [6]. Note 
also that given the overall mirror symmetry of our protocol, at 
most we only need £.Nt--,£-[{n+i)/2\ to be non-vanishing. 

To show this explicitly we must perform the direct calcula- 
tion, so first we introduce a convenient notation for our Clif- 
ford operators and their evolution. 

The operator A(f) = ^jX-^'jZ'j can be written as 

A(f) = 0'h3 with a = {xa^Za) = {xi,...,xn,Zi,...,Zn)- In this 
notation, Xr = Z^^t-^-xa-zb ^..^ ajjd the evolution is given 
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by 

T:A{t) — y A{t + l) 

a{t) — > a{t + l)=C.a{t) 

fit) — > f{t+i)^f{t)-{z{t)-m+m-L{m) 

where L{si,S2,.-,SnV = {s2,S3,...,Sn,0)''' , and 



C = 



r In 
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/O 1 ... 0\ 
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Vo 



10 1 
1 0/ 



A. The solution: 

Rather than try finding an elaborated numerical formula we 
shall follow the behavior for sample cases and from there infer 
the result, in fact the elements we depict are equivalent to take 
window ofthe...,Af/2±2,A^/2±l,A^/2stepsofaN element 
chain. Lets initially consider the even N case (e.g. = 8), 
with the notation —vx + uz^f. 



(7) 



It is straightforward then to compute, using e} = £, (—1,1) for 
all 

F'm-'F^[X{u)Z{v)],F'{T,)F' 

N+l 

= F^ (fg)-' [X{-u)Zi-v)], n P{£N+i-,n){m)F^ 

m=0 

^ ^-eS("[c"'S~+i-,»]-v,-"1c'"£«+,-,„].,) [x{u)Z{v)],. (8) 

This equation shows us that we can analyze the situation for a 
pulse at a single time m. We can ask ourselves what happens 
when we apply a pulse on time step + 1 — m: the phase 
factor in (8) 

^-v[C'"e„+i_„,].,-,+m[C'"e„+i_„,];, 

_ ^-l'[Ejv+l-ni{"l)]-V,-«[£/V+l-ra('«-l)].v, (-g-) 

_ ^-ve'^+l-"'(9(/+m-A'-l)-e(/-;n-l))-«E'^+'-"'(e(/+m-W-2)-e(/-m)) 

To get (10), we have used that a{t + I) = {x{t + l),z(r + 
1)) = (rx(f) +z{t),—x{t)) which with the boundary condition 
(x(0),z(0))==(-T,T) yields 

[x{t)]i^Q{l + t-N- l)-e(/-f - 1). (10) 

This means that a pulse in time N +1 ~m induces the trans 
formation 

f.N+1-m 
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X[u)Z{v) ^ "X(m)Z(v) for sites {;7z,A?+l-;7i}, 

^ lf^'-"'{-+>')x{u)Z{v) forsitesin [m+l,A?-i 

X(M)Z(y) otherwise. 

Our task is then to find a set of pulses which can generate 
actions on sites / and + 1 — / only. During our calculation we 
found it helpful to write a sample vector for every time step. 
Before proceeding we introduce some notation: we will call 
a /-plateau, M, (r) , a vector with equal non-vanishing compo- 
nents between row / and + 1 — ; only, additionally we will 
call an /-peak, Ri{r) as a vector with equal non- vanishing ele- 
ments on sites / and A^ + 1 — m only. Thus M,- (r) — M,+ 1 (r) = 
Ri{r), Mi[r) = M/(ar) and M/(r) +M/(r') = M,-(r + r'). So 
not to get the main argument lost in the subsequent section we 
want to restate what we seek: the problem is solved if we find 
a set of pulses reproducing an /-peak for any / as then Eq. (6) 
follows. Thus showing that universal Quantum computation 
is possible in a QCA fashion. We now proceed to show the 
solution. 



It is important to discriminate at this point between CV and 
the qudit case: when we have e5, then e can only take values 
within {0,(i- 1} for the qudit case and e M for CV. On CV 
our calculations are quite simple as we have such freedom that 
we can choose £ = 27t/(M + v) such that ^['=■^'1 = ^e(8;+Sjv+i-;) 
(turning 5/ into a /?/), however in the qudit case it is not al- 
ways possible. Nevertheless we can choose to apply pulses 
in different time-steps with different intensities (value of e) to 
get the desired result and thus in our example. 



Ra{u) 

Riiu^) 
Rx(u') 



= 84 = m(84 + 85) 

= mS's — (m + v)54 = M^(83 + Sg) 
= 11^82 — [u + v) {u8t, — v8n) ~ If' (82 + 87) 
= M^5i - [u + v){-i?-82 - V1183 + v^84) = i/(8i + 88) 

(11) 



now this wouldn't be successful if u" = d mod d for some 
s E {1,2, ...,A^/2}. However if such is the case then our anal- 
ysis can be simplified and choosing convenient intensities we 
would have, say u" = d. 
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Now the odd N case. 
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This case follows a similar recipe, 

7?4(v— m) =54 
^3(m(v — m)) = (v — m)5'3 — (m + v)5'4 
/?2(m^(v — m)) = m(v — m)52+ (v + m)(— (v — m)53 + v54) 
Ri (m^(v — m)) = ii^{v— u)S\ + (v + m)[— m(v — m)52 
+v(v — m)53 — v^54]. 

Now the critical case is when m'(v — m) = mod c/ for some 
s, however if that is the case then 

R2 iu'- ' v(v - m) ) = u'- 1 (v - ii)Si - m'" 1 (v - ii)Si+ 1 . 

as long as v — m 7^ mod d. When m = v mod d then we 
have the need of a global pulse ^o, which wasn't needed in 
previous cases to be able to generate a [^^yi]-peak, which is 
the main point to construct the rest of i-peaks, in our case we 
would have then that 

2Ra{u) = Sq + 2{S3+Si-S2) 

R^iu) = Sj,-2R4{u) 

Riiu) = 52-2/?4(m)-2/?3(m) 

Riiu) = Si-2{R4 + R3+R2) 

This shows that we can always find a set of pulses such that 

y (a, M, V, /) = n ■ e'«/2((l-C'^)^WZ(v) + (l-r'')(X(-»)Z(-v))t)A(;,/) 
^ -[-j''.g/a/2(l-cos(<|))){A:(«)Z(v)+Z(-v)A'(-M))A(;,/) 
g-tosm((t))/(Z(M)Z(v)-Z(-v)Z(-«))A(;,/) 

(12) 

with cos((|)) = so that all logical sites can be manip- 

ulated at will. Note that while some operations may not be 



executed on the middle q-site (which holds no information), 
e.g.M — v = modc/;M + v = mod t/, they can still be per- 
formed on the rest of the chain, i.e. the computationally valu- 
able q-sites. Here ^ and K are not completely freely chosen. 
Recall for CV, ^ = e' and for qudits ^ = e'-'^/''. The parameter 
K is a function of m and v determined by solving the equa- 
tions above, and is integer for qudits but can be real for CV. 
From (6), we would need [A^/2] + 1 — m non-vanishing pulses 
to construct an m-peak (/?„,), which leads to a total number 
4{N + 2) + —2m global operations. Solving the equations 
system for particular values of 11 and v simplifies the problem 
a lot and thus the number of global operations needed, e.g. for 
V = 0, Ri (m) = Si — Si+ 1 . 



B. Summary 

In the previous subsection we have shown that we are able 
to perform y(a,M,v, Z) and thus for whatever the value of 
k(m, v) we can run the protocol twice to either eliminate or re- 
inforce the Hamiltonian (X(m)Z(v) +Z(v)^X(m)^) noting that 
1 - - C"" e M and 5R(C^ - C,-"") = 0. This allow us then 
to apply any element of our operator basis on any logical site, 
thus achieving single qudit (qunat) arbitrary unitaries. To have 
complete universal quantum computation we only need now 
the ability to perform localizable entangling gates. 

To do so we can use a similar methodology to the one used 
in [6] for qubits. We can add |0) qudits (qunats) as ancil- 
las between logical sites, preserving the mirror symmetry of 
the encoding and thus preserving our previous results, such 
that a T-conjugated X;, TX+{s)iT-^ = ^(i),-i Z(-s),- (g) 
X{s)i+i +X\s)i-i iSiZ'^{-s)ii^X\s)i+i yields an effective 

logical e''^l'l®'^l'+'l'''^['1®^l'+i|' (or entangling gate (as |0) is sta- 
bilized by Z(i)). 

An alternate and more efficient way of implementing the 
entangling gate is the following: we perform the sequence to 
execute anX+(M) rotation on logical site one, e'"('^(")+'^(")^)i, 
then a time-shift 

pi, ^iaiX{ii)+X{uy^)i j^-m _^ gia(Z(M)„,®X(i()„,+i+Z(M),t,(8X(«)+^,) 

creates the required entangling gate between nearest neigh- 
bors. Note that this method does not need the extra |0) ancil- 
las, therefore reducing the overall number of gates needed to 
perform two q-site entangling gate. 

To show that this gate is enough for universal quantum com- 
putation we can just follow the argument in [3] for CV. For the 
qudit case we follow Ref. [ ], wherein the authors show that 
if one has a gate set consisting of arbitrary single qudit uni- 
taries and a gate which is diagonal in the computational basis, 
V \jk) = e'^j'' \ jk), which satisfies 

Qjk + Qpq T^^jq + ^pk, mod 271 for some y, A:, p,^, (13) 

then this suffices for exactly universal quantum computa- 
tion. Exact universality means that any unitary evolution 
can be obtained by a finite sequence of gates, as opposed 
to a dense of set of unitaries that cover the group. At our 
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disposal we have Uxx = (modulo sin- 

gle q-site unitaries), so we can have it's Fourier conjugated 

Uzz = e''^l'l®^l'+'l+4l®4+il' which has an action Uzz\jk) = 
^i2cos{^U+i^))\jk). Thus (13) turns into cos^{j + k) + 
cos ^{p + q) ^ cos ^{j + q)+ cos 2^ [p + k) mod 271, so for 
j = —k = 1 and —p = q = s we verify that 

2(1 +s)k 

2 ^ 2 cos fori^— 1, mode/. 

a 

Thus we have universal quantum computation. 

III. INITIALIZATION AND READOUT 

Readout of information in the chain can be done in an ar- 
chitecture independent manner using auxiliary states for each 
subsystem. First we consider the case of qudits {d finite). We 
employ two ancillary levels per subsystem, a level \a) that 
can be coherently coupled to at least one of the other d in- 
formation carrying states and another level |e) that couples 
to the state \a) but none of the information carrying states. 
The state |e) should couple to environmental degrees to al- 
low measurement by a classical readout and could represent, 
e.g. an optically excited state of an atom which decays emit- 
ting photons. High efficiency measurement of population in 
state \a) is possible if there is a closed cycling transition 
|a) <-*■ \e). To realize measurement on any qudit we simply 
adapt the formalism above but with all the Xj and Zj opera- 
tors acting on the d +1 dimensional subsystem spanned by 
{\0)j j , . . .\d — I) I , \a)j}. Then to measure population in 
the logical state \d — k)' of the qudit located at position / we 
apply the operator Xj', so that only the qudit located at site 
/ could have amplitude in an \a) state and then apply a uni- 
form pulse to couple \a) \e) and observe the presence of 
absence of a classical measurement output. By composition 
of these operations readout on any computational state can be 
performed. 

For measurement of CVs we could employ an additional 
degree of freedom per particle. Say we encode information in 
the X harmonic oscillator mode, and we have access to con- 
trol on a y harmonic oscillator mode initially prepared in the 
vacuum state for readout. The idea is to swap the state of the 
X mode (the data) of one CV into its y mode (the ancilla) and 
perform a global tomographic measurement on moments of 
the y modes of the CVs. This can be done in fact for CVs 
or for qudits using only global control and global measure- 
ments (see Fig. II). Our scenario is the following: (i) we use 



an ancilla chain initialized in the |0) state and (ii) demand that 
we can let the two chains (computational and ancilla) interact 
through a CZ gate (site ; of the main computational chain with 
site ; of the ancilla chain for all sites ;). With such conditions 
we can execute a swap gate between the q-site I of both chains, 
Swapd,,a, = {Fji®Fa,)CZdAFd,^®F,-^)CZd,aiF^^(^F^^) 
as shown in Fig. II. We have assumed that the ancillary CVs 
are prepared in the |0) state which are +1 eigenstates of the 
Z(s) operators. Such states are infinitely squeezed position 
eigenstates and are not physical (i.e. not normalizable). We 
can however achieve highly squeezed states by physically al- 
lowed global operations perhaps using the assistance of the 
associated data CV. Examples of such protocols include using 
using extremely short pulses of a standing wave potential [ 1 3] 
confining the spatial degree of freedom of the CV as demon- 
strated using trapped atoms in optical dipole potentials [ 1 4] 

We have considered here another chain, but we have pur- 
posely chosen only one of the chains to be capable of univer- 
sally controllable such that we can also think of every q-site 
having two degrees of freedom, e.g. an oscillator in y and x, 
but considering that q-sites can only be coupled though one of 
the degrees of freedom. 



IV. CONCLUSIONS 

In conclusion we have designed a protocol for univer- 
sal quantum computation with global operations for sub- 
systems encoding arbitrary finite or continuous variables. 
The data is stored in a mirror symmetric state over sub- 
systems encoding M q-sites aligned in one spatial dimen- 
sion. The overall requirements in our scheme are: (i) global 
nearest neighbors CZ gates, (ii) global Fourier pulses, (iii) 
global P{e„,) = X{—e)Z{E) pulses (iv) the set of Hamilto- 
nians {9l[X(fl)Z(fe)]}. For quantum computation with CV, 
condition (i) is met with the homogenous coupling H = 
81 1/ll' qi ® qi+u condition (ii) with H ^ g2Y!!^i{q] + pj); 
condition (iii) is met with H = Yli=i S'iiqi + Pi)^ ™d condition 
(iv) is met with Hamiltonians from the set {cosa)^,sina)^} 
with CO G M. We have also proposed a general scheme for 
readout using only global operations/measurements. 
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